preprint 



Elastic Tensor of Sr 2 Ru04 

Johnpierre Paglione, C. Lupicn, W. A. MacFarlane, J. M. Perz and Louis Taillefer 
Department of Physics, University of Toronto, Toronto, Ontario M5S 1A7, Canada 



O 

o 

X> 
PL, 



Z. Q. Mao and Y. Maeno 
Department of Physics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan 

(Dated: February 6, 2008) 

The six independent elastic constants of Sr2Ru04 were determined using resonant ultrasound 
spectroscopy on a high-quality single-crystal specimen. The constants are in excellent agreement 
with those obtained from pulse-echo experiments performed on a sample cut from the same ingot. A 
calculation of the Debye temperature using the measured constants agrees well with values obtained 
from both specific heat and Mossbauer measurements. 



c ■ 
o ■ 
o ■ 

5h ■ 

Or 

s : 

. 



i 

-a 
c 

o 
o 



> 

in 
o 

(N 
O 
(N 
O 
-i— > 

a 



i 

C 

O 

o 



PACS numbers: 62.20.Dc. 

Since the discovery of superconductivity in Sr2RuO4,0 
there has been much interest in the. possibility of novel 
pairing symmetry in this material. □ The availability of 
large, high-quality single crystals of Sr 2 Ru04 has allowed 
the full range of physical properties to be measured, with 
the prospect of identifying the order parameter symme- 
try and the mechanism of superconductivity. A thorough 
knowledge of the elastic properties is one important ele- 
ment in obtaining a complete understanding. 

The elastic constants of a material, which relate defor- 
mation to stress, are of interest because they are involved 
in fundamental solid-state phenomena: interatomic po- 
tentials, equations of state and phonon spectra. Further- 
more, thermodynamics links elasticity directly to quan- 
tities such as thermal expansivity, atomic, volume, Debye 
temperature, and Griineisen parameterB A determina- 
tion of the Debye temperature provides information on 
the phonon contribution to the low-temperature specific 
heat and the possible role of electron-phonon coupling in 
superconductivity. 

There are six independent second-order elastic con- 
stants Cij associated with the tetragonal crystal struc- 
ture of Sr2RuC>4. Using Voigt notation, they are ex- 
pressed as Cn, C33, C23, C12, C44, and Cqq. The con- 
stants for which i = j correspond to sound propagation 
in various principle crystal directions. When using a 
conventional time-of- flight, or pulse-echo, measurement 
technique, the determination of off-diagonal constants 
requires the troublesome measurement of sound propa- 
gation along non-principal directions, something which 
has not been done for Si'2Ru04. This difficulty does not 
arise when using resonance spectroscopy, which relies on 
a different technique to obtain the full elastic tensor of 
a material in a single measurement. In this article, we 
report the full elastic tensor of Sr2Ru04 obtained by res- 
onance spectroscopy. 

The single-crystal sample of Sr2Ru04 was grown by 
the traveling solvent floating zone techniques and found 
to have a superconducting transition temperature of 1.37 
K as determined from magnetic susceptibility, which is a 
good indication of its high purity and quality. The sam- 
ple was aligned using Laue x-ray diffraction, cut into a 



rectangular parallelepiped using spark erosion, and pol- 
ished to dimensions of 2.88(1) x 1.99(1) x 0.86(1) mm 3 . 
The sample mass was 0.0292 g, yielding a density of 
5.92(8) g/cm 3 , which agrees with the previously reported 
value of 5.918 g/cm 3 . □ 

Measurements were made at room temperature us- 
ing resonant ultrasound spectroscopy (RUS)Ja This tech- 
nique abandons the plane-wave approximation used in 
pulse-echo experiments and instead uses the normal 
modes of vibration of a solid specimen of known geome- 
try, crystal symmetry and density to deduce the complete 
elastic tensor in a single measurement. Thus, it allows 
for relatively easy determination of off-diagonal elastic 
constants, which are generally impractical to measure in 
non-cubic symmetries using the pulse-echo technique. 

The acoustic resonance frequencies of a single-crystal 
specimen can be calculated given the dimensions, mass, 
and elastic constants. The key to the RUS technique lies 
in the ability to determine unknown parameters (in this 
case the elastic constants) from knowledge of the reso- 
nance frequencies, which are readily determined experi- 
mentally. There exists no analytic method of perform- 
ing this calculation, and so the unknown parameters are 
determined via a computational fitting procedure which 
uses an iterative algorithm to match resonance frequen- 
cies calculated analytically with those measured exper- 
imentally. This calculation requires the input of mea- 
sured dimensions, mass, and an estimated set of elastic 
constants to start the fitting iteration: the elastic con- 
stants are left as adjustable parameters to be determined 
by minimization of error between the measured and cal- 
culated frequencies. Because the largest source of error 
in this experiment was the determination of the sam- 
ple dimensions, they were also left as free parameters in 
the fitting routine, allowing for a comparison between 
measured dimensions and those determined from the fit. 
As a test, the RUS apparatus and method were used to 
obtain the complete elastic tensor of the heavy-fermion 
material UPt3 (which has a hexagonal crystal structure). 
This gave the values C l2 = 1.44(2), C 23 = 1-70(2), 
C 33 = 2.93(2), C44 = 0.3794(1 jk and C 66 = 0.83 6 7(9) 
(all in units of 10 12 dynes/cm 2 ) p which are in excellent 
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FIG. 1: A portion of the resonance spectrum of a rectangular 
parallelepiped of S^RuCU obtained using the RUS technique. 
The inset shows a schematic of the RUS apparatus. 



agreement with those measured by de Visser et al. using 
the pulse-echo technique £j 

The resonance frequencies of the Sr2Ru04 sample were 
determined at room temperature using the apparatus 
shown schematically in the inset of Fig. |l|. By holding 
the sample between two piezoelectric transducers, where 
one transducer is used to drive the sample with mechani- 
cal vibrations and the other to detect the mechanical re- 
sponse of the sample, a frequency sweep allowed a rapid 
measurement of the resonance modes. It should be noted 
that the only assumption made in the RUS technique lies 
in the ability to observe the natural resonances of a solid 
body under free boundary conditions, allowing a proper 
comparison between measurement and calculation. The 
rectangular parallelepiped sample was held only lightly 
by its corners to minimize suppression of any resonance 
modes, thus making the free boundary approximation 
appropriate. 

Measurements up to ~ 3 MHz provided more than 
sixty measured resonance frequencies. The central por- 
tion of this frequency spectrum is shown in Fig. The 
first thirty measured frequencies are listed in Table | 
along with those calculated analytically using a program 
developed specifically for the rectangular parallelepiped 
sample geometry.Q Note that all of the first thirty cal- 
culated resonance frequencies were detected, and only 
three of the next thirty were undetectable. The error is 
listed as the percentage difference between measured and 
calculated values. It should be noted that in RUS mea- 
surements it is commonly found that the first one or two 
modes typically have a much larger discrepancy between 
measured and calculated values than the average, as is 
seen for the first mode in Table B, and are therefore not 
weighted in the fit. The dependence of each mode on the 



elastic constants was also estimated using this program 
in order to indicate the pure (shear/compressional) or 
composite nature of each resonance mode, and are listed 
as the normalized values df / dCij in Table |. For exam- 
ple, the n — 5 resonance mode can be seen to depend 
mainly on the constant Cqq , and is thus primarily a pure 
shear mode. Identification of pure resonance modes can 
be of use in tracking the behavior of individual elastic 
constants as a function of temperature or magnetic field. 

The elastic constants used as initial fitting parameters 
were estimated using a combination of available values for 
S^RuOjiini-J and the isostructural material La2CuOji3, 
and subsequently adjusted through numerous fitting it- 
erations to obtain the best fit. The quality of each fit 
was calculated using the root-mean-square error orms 
between measured and calculated resonant frequencies, 
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where n is the mode number and N is the total number of 
modes weighted in the fit. The fitting procedure resulted 
in a ctrms value of 0.19% between the measured and cal- 
culated frequencies; such a small value is an indication 
of an excellent fitfl The fitted dimensions agreed with 
the measured values to within error, also confirming the 
validity of the fit. The elastic constants calculated using 
this fit are compared in Table || with values obtained for 
Si'2Ru04 from various pulse-echo experiments and with 
those obtained for La2CuC>4 using the RUS technique. 
The absolute accuracy of the constants obtained in this 
study was based on either a quality-of-fit estimate, ob- 
tained from a measure of the stability of the minimiza- 
tion calculation used in the fitting routine, or on errors 
in density measurement, taking the larger of the two as 
the conservative approximation. Specifically, the errors 
on Cn, C33, C23, and C12 were obtained from quality- 
of-fit estimates, and those on C44 and Cqq were based 
on the variation in calculated dj values obtained from 
fits using upcer and lower dimensional bounds as input 
parametersJlj 

The elastic constants reported by Lupien et al. were 
measured at low temperature (< 4 K) using the pulse- 
echo technique, on samples cut from the same ingot used 
in this studyJlil A typical pulse-echo experimental appa- 
ratus is shown schematically in the inset of Fig. ^. With 
this technique, sound velocities are calculated using the 
length between two parallel faces of an oriented crystal 
and the measured time delay between echoes of the initial 
pulse, such as those shown in Fig. ||. By solving the usual 
Christoffel equations, the elastic constants are obtained 
from sound velocities measured for various propagation 
and polarization directions and the mass density.tj For 
example, the shear elastic constant, 

Cm = P«[ioo/oio]) ( 2 ) 

is calculated from the mass density p and sound velocity 
U[ioo/oio] (propagating in the [100] direction and polar- 
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TABLE I: Comparison between measured and calculated values of the first 30 resonance frequencies of the S^RuCU sample. 
Included are the normalized df/ddj values, which indicate the nature of each resonance mode n. 
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"Calculated for comparison using measured values of Cn and C12 . 
b Cij values extrapolated to 300 K using data from Ref. ha. 
c Measured at 300 K. 



ized in the [0101 direction). All of values from Ref. |l4| 
shown in Table |ll|, except for C12, are proportional to the 
square of various velocities directly measured in differ- 
ent propagation and polarization directions (the value for 
C12 is given as an average of two values calculated using 
two different sound velocities and the value of Cn). By 
extrapolating these low-temperature values to 3JX) K us- 
ing temperature dependencies found elsewhere, 113 agree- 
ment is found to within experimental error for C44 and 



C66, and to within 10% for Cn, Cn and C33 when com- 
pared to the values in this study. The various elastic 
constants reported by Matsui et al. (Ref. 11, |l5|) and 
Okuda et al. (Ref. |l6|) are a factor of ~ 2 less than these 
values, which is much greater than the ~ 6% change in 
the constants observed from low to room temperature. 



The Debye temperature 9jy is estimated from the elas- 
tic constants by using the following expression for a 



4 



6 




50 100 

Time (us) 



FIG. 2: Echo pattern for a transverse sound mode corre- 
sponding to Cm, with propagation in the [100] crystal axis 
direction and polarization in the [010] direction of S^RuCU, 
obtained using the pulse-echo technique on a sample of length 
3.98 mm (from Ref. 0). The inset shows a schematic of a 
typical pulse-echo setup. 
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where h and ks are Planck's and Boltzmann's constants, 
respectively, N/V is the atomic density, and J is the 
series expansion of an integration over sound velocities 
in all crystal directions, written in terms of the elas- 
tic constants (see Eq. 32 in Ref. [l?]). Using Eq. 3, the 
full set of elastic constants measured in this study gives 
9d = 465(5) K (this value increases by ~ 5 K when 
using values extrapolated to low temperature in a man- 
ner similar to that used in Table [n] for Ref. [l4j). This 
elastic measure of 9d compares well with the value of 
410(50) K extracted from the T 3 phonon contribution to 
specific heatjia and the value of 427(50) K obtained in 
Mossbauer measurements,!^ where 6d is extracted from 
the temperature dependence of the Debye- Waller factor. 

In conclusion, the complete set of elastic constants was 
obtained for the tetragonal crystal Sr2Ru04 using the 
measured resonance frequencies of a small rectagular par- 
allelepiped sample. Excellent agreement was found with 
values obtained from pulse-echo experiments performed 
on a sample cut from the same ingot. The Debye tem- 
perature determined from the measured constants was 
found to compare well with values obtained from other 
measurements. 
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